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Abstract 

We test the positivity property of the chiral perturbation theory 
(ChPT) pion-pion scattering amplitudes within the Mandelstam trian- 
gle. In the one-loop approximation, 0(p 4 ), the positivity constrains 
only the coefficients 63 and 64, namely one obtains that 64 and the 
linear combination 63 + 364 are positive quantities. The two-loops 
approximation gives inequalities involving all the six arbitrary param- 
eters entering ChPT amplitude, but the corrections to the one-loop 
approximation results are small. ChPT amplitudes pass unexpect- 
edly well all the positivity tests giving strong support to the idea that 
ChPT is the good theory of the low energy pion-pion scattering. 
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Key-words: Pion-Pion Scattering, Chiral Perturbation Theory, Pos- 
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1 Introduction 

Chiral perturbation theory (ChPT) is considered a low-energy effective ap- 
proximation of QCD, in particular it provides a representation of the elastic 
pion-pion scattering amplitudes that is crossing symmetric and has good 
analyticity properties. In a seminal paper [0] Gasser and Leutwyler devel- 
oped ChPT which allows one to compute many Green functions involving 
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low-energy pions. It is well known that the physical pion-pion scattering am- 
plitudes can be expressed in terms of a single function A(s,t,u) whose form 
was obtained as a series expansion in powers of the external momenta and of 
the light quark masses. The first term of the series was given by Weinberg 
[§], the second by Gasser and Leutwyler [L] and only recently a two- loop 
calculation 0, f|] was obtained. In this approximation the function A(s, t, u) 
has the following form 



A(s, t, u) = a(s — 1) + a 2 bi + b 2 s + b 3 s 2 + b±(t — u 



F^(s) + G {1 \s,t) + G w (s,u) 



+ a° 



b 5 s 3 + b 6 s(t — u) 



+ 



+ 



'F®{s)+GV>(s,t) + G®\ 



s,u 



0[a A 



where a = (M^/F^) 2 , M n is the mass of the physical pion, F n the pion decay 
constant, s, t, u are the usual Mandelstam variables, expressed in units of the 
physical pion mass squared M 2 

s = ( P i+ P2Y/MI f = (pi - P2Y/MI u = (pi - PsY/M 2 , 

F^(s) and G^\s,t) are known functions and 6j, % — 1, . . . , 6 are arbitrary 
parameters which cannot be determined by ChPT JL], |3], |]]. In any realistic 
comparison with experiment we have to provide some numerical values for 
all these parameters obtained from other sources. One hopes that by using 
unitarity this can be done although, until now, no program for implementing 
this property has been presented. The common belief is that imposing uni- 
tarity is not a simple matter since its implementation in one channel destroys 
crossing symmetry in other channels. However there is a weak form of uni- 
tarity, the positivity of the absorptive parts, which is a linear property and 
which can be imposed. This property was used thirty years ago to obtain 
constraints on the 7r°7r° s-wave partial amplitude, fo{s), in the unphysical re- 
gion < s < 4 and on the (i-wave scattering lengths. These constraints were 
useful because at that time almost nothing was known about the explicit form 
of the scattering amplitudes and they were used in testing models for pion- 
pion partial-wave amplitudes. The advantage of ChPT is that it furnishes an 
explicit form for the pion-pion scattering amplitudes whose unknown part is 
contained in a few numerical coefficients. Thus is of certain interest to see 
how these properties reflect into constraints on the bi coefficients entering 
Eq.(l.l). 
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Beginning with the paper B, Martin has used the positivity, analyticity 
and crossing symmetry to obtain constraints on the 7r°7r° s-wave partial am- 
plitude, fo(s), in the unphysical region < s < 4; a few of them have the 
following form || 

/ (4) > / (0) > / (3.15) / (0) > \ £ f (s) ds 

F(0,2(l + i=))>/ (2(l + i=)) (1.2) 

where F(s, t) denotes the 7r°7r° elastic scattering amplitude. A more complete 
set is found in ref. [[7]. The most elaborate form of these constraints is the 
following result: the 7r°7r° s-wave, fo{s), has a minimum located between 
1.218989 < s < 1.696587 §, 0, |, |, [H| 0, Q and this result can be 
improved only by unitarity. These results can be translated into constraints 
on the parameters bi entering ChPT pion-pion scattering amplitudes. As 
we will see later all the above inequalities are equivalent in the one-loop 
approximation, 0(p 4 ), to a single constraint on the coefficients bi whose 
typical form is 

37 

63 + 36 4 > 



1920vr 2 

the only difference being the numerical value appearing on the right hand 
side. This is the explanation of the inefficacy of these apparently distinct 
constraints which was observed from the beginning by people constructing 
models for pion-pion partial-waves. 

In this paper we work only in the unphysical region \s, t, u\ < 4, i.e. there 
where the amplitude (1.1) is considered to be a very good approximation 
to the true amplitude. Other approaches use information from the physical 
region to obtain constraints on the same parameters bi ]13, [14], [H|. In the 



one-loop approximation, 0(p A ), the positivity property constrains only the 
coefficients 63 and 64. By taking into account the 0(p 6 ) contributions one 
gets constraints involving all the six parameters entering Eq.(l.l). 

By using the unitarity bounds on 7r°7r° scattering amplitude in the un- 
physical region |lj| one get upper and lower bounds on some linear combina- 



tions of the parameters 6j. These unitarity bounds are not very constraining; 
to see this we give that one obtained from the bound on F(2,0), where as 
above F(s, t) denotes the 7r°7r° amplitude. The bound is —3.5 < F(2, 0) < 2.9 
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and it is equivalent to the following lower and upper bounds 



-3.5 x 32tt < a + a 2 ( 3&i + 46 2 + 86 3 + 864 + — ^ ) + a 6 [166 5 + 166 6 + 

327T^ 



(4-tt) /56i 62 H63 56 4 \ 965 251 41 



7T 2 \ 16 2 12 12 / 3456vr 4 3456tt 3 6144vr 2 



< 2.9x32tt 



Because of the factor 327r appearing on left and right hand side the bounds 
are not very strong and we will not consider them here. 

The physical isospin amplitudes F 1 can be expressed in terms of the single 
function A(s, t, u) as follows 

F°(s,t,u) = 3A(s,t,u)+A(t,u,s)+A(u,s,t) 
F 1 (s,t,u) = A(t,u,s) - A(u,s,t) 
F 2 (s,t,u) = A(t,u,s) + A(u,s,t) 

where A(s,t,u) is given by Eq.(l.l). 

Having only three independent amplitudes one gets only three indepen- 
dent constraints since the crossing symmetry is an exact symmetry for the 
ChPT amplitudes. The construction of our positivity constraints is outlined 
in the next Section where we present an over-determined system of con- 
straints. Their implications on the coefficients 6j are discussed in Section III. 
The paper ends with Conclusion. 



2 Positivity constraints 

Let F 7 (s,t) denote the ittt scattering amplitude with isotopic spin / in the 
s channel. In matrix notation F(s, t) satisfies the following crossing relation 

[0 

F(s,t) = C st F(t,s)=C su F(u,t) 
where the notations are 

F(s,t)= F 1 (s, t) 

\F\s,t) ) 
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/ 1/3 1 5/3 \ / 1/3 -1 5/3 \ 

C st = 1/3 1/2 -5/6 , C su = -1/3 1/2 5/6 
V 1/3 -1/2 1/6 / V 1/3 1/2 1/6 / 

From the results of axiomatic field theory we know that the amplitudes 
F J (s,t) satisfy fixed-t dispersion relations with two subtractions |18[ for 
\t\ < 4. We may write them as 



F(a, t) = C st [a(t) + (a - u)b(t)] + - / — + A(s, t) 

(2.1) 

where A(x,t) is the absorptive part of F(s,t) and the subtraction constants 
are of the form 



a(t) 




b(f) 




due to crossing symmetry. 

In the following we shall consider that s, t, u take values in the unphysical 
region \s,t,u\ < 4. We calculate the difference 

F(a,t)-F(s a ,t) 

and we are looking for those combinations of isospin amplitudes for which 
this difference does not depend on the subtraction constants. From Eq.(2.1) 
we find 

(F(s, t) - F( Sl , t)) = 2C st b(t) + f(A) (2.2) 

s — Si 

where f(A) denotes the complicated term containing the integration over the 
absorptive parts. The first term on the right hand side of Eq.(2.2) is 



c st h{t) 



( b(t) 
\b{t) 

V -W) 



The last relation shows that there are three combinations of isospin am- 
plitudes for which the difference (2.2) have no dependence on the subtraction 
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constants. They are F° + 2F 2 , F l + F 2 and F° — 2F 1 and we shall denote 
them as Fj, % = 1,2,3 in this order. The first one is the well known 7r°7r° 
elastic amplitude. One easily obtains from Eq.(2.2) the relation 

s-si)(s-wi) f°° (2x + t - 4) Ai(x, t) dx 



F,(., f )-fi(.„t)= "— ^-"V 

71 J A 



A (x — s)(x — Si)(x — U)(x — Ui) 



1, 2, 3. From this relation we get 

dFi(s,t) s-u f°° (2x + t-4)Ai(x,t)dx 



d-S TT J A (x — s) 2 (x — U) 2 

Because the absorptive parts A\ and A 2 are positive we find that 

1 dFi(s,t) 



u ds 



> i = l,2 



The third combination involves the absorptive parts A°(x,t) — 2A 1 (x,t) 
whose sign is not defined and we cannot say anything about the sign of 
the derivatives of F^(s,t). The precedent relations show us that on the line 
s = u, Fi(s, t), i — 1,2 attain their minimum values. Indeed we obtain from 
them the second derivatives 



d 2 FAs,t) 2 roo ( 1 1 \ 

;V - = - / 7 + 7 v )Mx,t)dx>0, < = 1,2 2.3 



which are positive definite, implying that the functions Fi(s,t) have a mini- 
mum on the line s — u. From the last relation we obtain also 



9 2n_1 F i (s,*) (2n-l)\ ( 1 



5s 2n-l 



7T 



1 



(x — s) 2ri (x — u) 



2n 



Ai(x, t) dx = 



(2n-l)\(s-u) [™ 



roo 

y 4 [(x - s ) n + (x - u ) n ] 



(x - s) n + (x — s) n (x - u) + . . . (x - u) 
(x — s) 2n (x — u) 2n 

In this way we obtain the set of positivity constraints 



n-l 



Ai(x, t) dx 



i a fc -Ji(«.0> 0aaJ g fc f.(«.0> (2 . 4) 



S — M 



1 — 
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i = 1,2, n = 1,2,... 

A first remark is the following, if the positivity constraints have to be 
fulfilled it is sufficient to test them only on the line s — u, i.e. 2s + 1 — 4 = 0, 
where the functions Fi(s,t) attain their minimum values. In this way we 
have only one free parameter, < \s\ < 4, and on this line the odd and even 
derivatives give the same information. From the point of view of computation 
it is simpler to work with even derivatives. 

Up to now we have obtained two constraints given by Eq.(2.3). Because 
we have three independent isospin amplitudes it follows that we can obtain 
another one at most. 

The positivity constraints can be imposed even on the isospin amplitudes 
themselves. This can be easily seen from the relation (2.1) for F 2 (s, t) which 
after derivation gives 



d 2 F 2 (s,t) 2 , . 

'ax 



ds 2 



2 r 

7T J 4 



(x — s) 3 6 (x — u) 3 



' ' A 2 {x,t) + 



' 1 l^M + lj-^A^t) 



3(x-u) 3 v ' ' 2(x-u) 3 ' 

The right hand side of the previous relation is a positive quantity and by 
iteration we obtain that 

„=!,,,... (,5> 

Unfortunately the numerical calculations show that this relation is not in- 
dependent of the previous two ones. An other way to obtain them is to make 
use of the Gribov-Froissart representation for the partial wave amplitudes. 
One writes dispersion relations for the isospin amplitudes, the subtraction 
constants being given by the s- and p-wave partial amplitudes, and one finds 

1 r°° 

F T (s, t) = f*(s) + - / A f (x, s)g(x, s, t) dx 1 = 0,2 (2.6) 

7T JO 

where 

g(x, s, t) = + + ln(l + 

x — t x — u 4 — s x 

For I — 1 we can write similarly 

Of 1 roo 

F 1 (s, t) — 3/ 1 1 (s)(l + t) + -/ A\x,s)h(x,s,t)dx (2.7) 

s — 4 TC Jo 
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where 
h(x, s, t) 



x — t x — u 



1 6(2* + s - 4) 



(4 - S )s 



2x + s-4, s-4 x 
, I + — — ) + 2 



The absorptive parts entering (2.6)-(2.7) are the t-channel ones and in the 
following we make use of the t < — > u crossing symmetry. From the relation 
(2.6) we find analogous formulas to (2.4), namely 

v ' ; > and ^ ' ; > (2.8) 



t — u 



Q t 2n- 



1 — 



dt 2r - 



n = 1,2, . . . , 7 = 0, 2, which are not numerically independent of the previous 
ones. More interesting is the relation (2.7) which can be written as 



t — u 



3/M 



-/ A\x,. 

7T JO 



= F*(M) = ^j + 

S — 4 vi ju 

2x + s- 4, s - 4 



(x — t)(rr — w) 



(4 



4 



-Zra(l + 



x 



)+2) 



dx 



This relation provides us with another independent relation. Because F 1 (s, t) 

is antisymmetric in t < ► u, F 1 (s,t) is an analytic function within the 

Mandelstam triangle. Deriving it with respect of t we obtain 



dF^t) t-u 



A 1 (x, s) 



dx 



dt it Jo (x — t) 2 (x — u) 2 

From this relation we obtain similar formulas to Eq.(2.4), namely 



t — u 



Q t 2n- 



Q t 2n 



(2.9) 



Finally the positivity may also be expressed as the positivity of the partial 
wave amplitudes 



f!(s) 



4 - s h 



A*(x, s)Qi 



2x 
4- s 



l)dx 



(2.10) 



for / > 2 inside the unphysical region < s < 4, but the numerical calcula- 
tions show that these last constraints are weaker than those derived above. 
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3 Numerical Results 



In the previous section we derived a complete set of positivity constraints. 
In an exact theory many of them are consequence of the others as we will see 
later. Because ChPT does not completely specify the amplitude and on the 
other hand the power of different constraints is not the same it is useful to 
derive as many constraints as possible. Since all of them have to be satisfied 
we will select the strongest one in every case. 

To test the method we worked first in the one-loop approximation, 0(p A ), 
i.e. we retained terms up to a 2 in Eq.(l.l). In this order the obtained con- 
straints do not depend on the value taken by a; in the two-loop approximation 
the constraints will be linear in a. We consider first the constraints on the 
7r°7r° scattering amplitude. 

The first two constraints (1.2) on the 7r°7r°, s-wave /o(4) > /o(0) > 
/o(3.15) are equivalent to 

9 29 

h + 3b 4 > -1.64 xlO- 2 and b 3 + 3b 4 > -1.47 x 1(T 2 
lUz4 oo47T z 

respectively. The third relation /o(0) > \ f£ fo(s) ds is equivalent to 

h + 36 4 > -4.64 x 1(T 4 

the strongest result being the last one. Already the last relation (1.2) fur- 
nishes a better result, the above combination of coefficients gets positive 

b 3 + 36 4 > 7.36 x 10" 4 
The derivative of the s-wave has the form 

fo(s) = l(h + 3b 4 )(5s-8) + h(s) 

where the function h(s) has a long expression which we do not write it here. 
Since 5s — 8 < 0, for s < 1.6, the upper and lower bounds on the derivative 
/g(s) describing the position of the minimum are equivalent to a single lower 
bound on the combination b 3 + 36 4 . The best result is obtained at s = 
1.696587 and is 

b 3 + 3b 4 > 6.67 x 10 -3 
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The last result is the strongest constraint upon the combination 63 + 364 
obtained from inequalities satisfied by the 7r°7r° s-wave within the unphysical 
region < s < 4. We have given all the above results to understand why 
these constraints were easily satisfied by the phenomenological models for 
partial-wave amplitudes constructed in the past years; satisfying a few of 
them the others are automatically fulfilled. 

A similar analysis with analogous results was done in ref. 
the game the s- ans d-waves. 

Stronger constraints are obtained from the positivity of the second deriva- 
tive of the full amplitudes. From Eq.(2.3) we find for % = 1, i.e. the 7r°7r° 
amplitude again, a relation of the form 

6 3 + 3& 4 + /ii(s) > 

where hi(s) is a decreasing function for s < 4. Thus a problem arises, at 
what point has to be considered the above relation. We decided to limit the 
range of s within the interval \s\ < 4 since we are aware that the amplitude 
(1.1) is only an approximation of the true amplitude, approximation truly 
not valid for values s > 16 in the physical region, s = —4 is equivalent to 
t = 12 in the physical region of the t-channel. One gets 

37 

h + 36 4 > o ~ 1-92 x 1(T 3 , for s = and 

3 ~ 1920vr 2 

63 + 36 4 > 8.28 x 1(T 3 for s = -4 

which is stronger than the previous inequality. In the following we will list 
numerical values only at s = 0, the numerical values at s = —4 being not 
very different although they are a little better such as the previous relations 
show. 

We will work now in the two-loops approximation, 0(p 6 ), and consider 
only the constraints derived in Section 2 which give the strongest results. 
The constraints have the form 

6 

c-(s,a) bi + f k (s,a) > 

i=l 

where k = 1,2,3 labels the independent constraints. As we said before 
there are only three possible constraints from the positivity of the absorptive 



including in 
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parts, because we have only three independent amplitudes and the crossing 
symmetry is an exact symmetry for the ChPT amplitudes. In the previous 
section we derived six constraints but only three are numerically independent. 
The above inequalities are obtained for every value of s within the unphysical 
region < \s\ < 4. Each inequality defines a half space where can live the 
parameters fej. Thus the true result would be that obtained by constructing 
the intersection of these half spaces. Unfortunately the envelope cannot be 
found in an analytic form the functions fk(s,a) being very complicated. 
The constraint (2.3) for i — 1 and s = is equivalent to 



b 3 + 3b 4 



37 



1920tt 2 



+ a 



7h 



2b, 



+ 



320tt 2 60tt 2 45vr 2 180tt 2 



+ 16& 6 



367 



+ 



6869 



Making a 
For i = 



552960tt 2 1658880tt 4 
-> one gets the one-loop result. 
2 we find the relation 



> 



31 



5760tt 2 



bi 



436 2 h 236 4 



(3.1) 



240tt 2 2880tt 2 24tt 2 180tt 2 



67 



707 



276480tt 2 331776tt 4 
The inequality (2.5) for n = 2 gives 



> 



46 fi 



(3.2) 



b 3 + 5 6 4 - 



173 



5760tt 2 



+ a 



13 61 



67 6 2 



23 6, 



960tt 2 1440tt 2 180tt 2 4tt 2 



+ 24 6 6 - 



11 



+ 



13939 



> 



61440tt 2 1658880tt 4 „ 
and it is easily seen that it is a consequence of the previous two ones being 
the sum of the first and of the second one multiplied by 2. The relations 
(2.6)-(2.8) give us, in principle, three new inequalities, but only one will be 
independent of the first two already obtained. For I = we obtain 



63 + 764 



47 



11527T 2 



+ a 



11 6 2 19 63 



91 64 

1927T 2 ~ 144tt 2 ~ 907 2 ~ 1807T 2 



+ 32 6 6 + 



11 



169 7003 
552960vr 2 + 550960vr 4 



> 



which again is a linear combination of the first two ones. 
For 1 = 1 the relation has the form 



11 

+ a 



2688vr 



h 17b 2 15 16 3 6 536 4 

+ TTTTT^ - 7777777^7 " TTTTTtA + h + h + 



224tt 2 1344tt 2 2016tt 2 2016tt 2 



37 4111 
+ 



> (3.3) 
215040tt 2 290304vr 4 J ~ v ; 

This is the third independent relation as can easily be seen because the 
one-loop approximation gives a positive number independent of 6j. More 
important is the fact that the function f 3 (t, a) appearing in the relation (3.1) 
in the one loop approximation is positive over (presumably) all the negative 
real axis which proves that the positivity is very well satisfied even by the 
lowest approximation of the ChPT amplitudes! 

For 1 = 2 the inequality is 



. 7 49 

03 + Oa + a 

5760tt 2 



29b, , 19b 2 7b 3 A7b A 

+ 77777 n + 7777; n + 8 6 " 



9607T 2 14407T 2 1807T 2 1807T 

127 67 



, > (3.4) 
110592vr 2 552960vr 4 J ~ v ' 

We have written all the inequalities since they were useful in checking the 
calculations. 

A first remark is the following, the coefficient 65 appears only in the 
inequality (3.3). Thus we can say that it is unimportant and make it vanish 
also in the amplitudes! It is true that the above inequalities have been 
obtained by an extremal property, they are calculated on the line where the 
corresponding amplitudes are taking their minimal values. This may be a 
suggestion that the physical partial waves satisfy also an extremal princ iple 
which has to be found. This is also supported by the findings of Wanders, 



who could not obtain a reliable value for this parameter W5\. What the above 
results suggest is that a good determination of 65 can be obtained only from 
1=1 data. 

We have tested how the parameters bi found in literature compare with 
the inequalities. Unfortunately there are only two papers that gives values for 
all bi pm . The values given by Bijnens et al. are in the domain allowed by 
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Eqs.(3.1 )-(3.4), but the values obtained by Ananthanaryan strongly violate 
the inequality (3.4), the previous ones being satisfied. The values obtained by 
Knecht et a/.|Q for the last four parameters can be used to obtain constraints 
on the bi an d 62 but the allowed domain is rather large; the same is true for 
the values given by Wanders jnj. 

Using Roy equation analysis of the available tttt phase shift data Anan- 
thanaryan and Buttiker[21] obtained values for the chiral coupling constants 
Ii, I2 but their results are not easily translated into constrain ts on the bi 
coefficients We have tested also the inequalities (2.4), (2.8), (2.9) and (2.10) 
for a few values n > 3 and I > 2 respectively and we found that they are 
very well fulfilled. 



4 Conclusion 

We have tested the positivity properties of the ChPT pion-pion amplitudes 
and we have obtained a number of inequalities which express this property. 
We conclude that the pion-pion amplitudes given by the relation (1.1) satisfy 
this property unexpectedly well. In the 0(p ) approximation the positivity 
implies two constraints: 64 is a positive quantity and so is the combination 
b 3 + 36 4 . As concerns the positivity of the 1 = 1 amplitude this was tested 
numerically up to t = — 10 5 where the derivative is still positive. Including 
0(p 6 ) contributions one get constraints involving all the six parameters bi, 
but the corrections to the 0{p i ) results are small which support the idea 
that the expansion (1.1) is the best candidate for the true amplitude. It 
might seem surprising but we consider that the main result of this study 
is the conclusion that the most important parameters to be determined are 
63 and 64. This is a consequence of the good properties near threshold of 
the Weinberg approximation together with the very powerful property of 
positivity of the scattering amplitudes. Let us remind that this property 
was essential in deriving the analyticity domain of pion-pion amplitudes 



This means that the amplitude (1.1) in which all but 63, 64 are zero will give 
a fair description of the low energy phenomenology up to about 6-700 MeV 
and the contribution of the other coefficients will be seen at higher energies. 
In conclusion a good determination of the above two parameters will be a 
good starting point in the comparison of the theory with experiment. Work 
on this line is in progress. 
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